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5_i I Abstract: 



For an accurate description of the polarized deep inelastic scattering at low x including the logarithmic correc- 
tions, ln^(l/2;), is required. These corrections resummed strongly influence the behaviour of the spin structure 



Ph. 
< 

^ functions and their moments. Results of the work of J. Kwiecinski and myself on this problem are reviewed, 

> , 

(yQ , 1. Introduction 

t:;;^- ' The problem of identifying the spin components of the polarized nucleon has been at- 

I tracting the attention of the high-energy community since many years The data 

I obtained in 1988 by the EMC collaboration [3] showed that the total participation of quarks 

' in the proton spin was very small. This contradicted theoretical predictions, obtained with 

Q-i' 

I ' the well-founded Ellis- Jaffe sum rule That sum rule connected the moments of quark 

Qh' 

Q I distributions to the nucleon axial coupling constants. Following that rule, the quarks should 
(— ! I 

' participate in about three-fifth of the total nucleon spin with the parton quark model. The 
discrepancy between the theoretical expectations and the experimental data for the polarized 

X : proton ha. been often referred a. the "puzzle of the proton ,pin" 0. 

c3 

Since 1988 several observations contributing to this problem have been made. New ex- 
perimental data were collected. They showed that neglecting the participation of the strange 
quarks at the Ellis- Jaffe sum rule was not correct |31I2] • 

A significant progress in the theoretical analysis of the nucleon spin structure function 
has been also made. A derivation of gi in the subleading approximation of the quantum 
chromodynamics (NLO DGLAP) showed a difficulty with the unique defining the quark 



component of the polarized nucleon. In the standard MS factorization scheme at the 
NLO QCD |ni|7j the first moment of the singlet parton distribution, AT,{Q'^), is not con- 
served. The singlet moments calculated at different factorization schemes differ by a term, 

as{Q'^) / dxAg{x,Q'^) ~ 0{a^), which can be large even at oo. This is a direct con- 

Jo 

sequence of the axial anomaly [HlISllS] . Therefore the differences between the quark moments 
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at different scfiemes can be large. Tfiis makes tlie pliysical interpretation of tliese moments 
difficult. 

It is also not clear how the orbital angular momentum, L, participates in the total nucleon 
spin P]. Again, there is a problem with the unique definition of this momentum. The total 
angular momenta of quarks and gluons, Jq and Jg, are well defined and gauge invariant. 
However, separation of their spin and orbital components is not unique. It depends on the 
gauge transformation used. 

New experiments emerge, and they will possibly help explaining the puzzle of the nucleon 
spin. The data from the region of low values of Bjorken x, x < 10~^, will be of special 
importance. Theoretical predictions show that at low x the structure function, gi{x,Q'^), 
is influenced by large logarithmic corrections, /n^(l/x), jlUllllj . As a consequence, large 
contributions to the moments of the structure functions from this region are expected jlUl 
Illlll2lll51ll4lll5lll6j . Including the region of low x into the experimental analysis will improve 
the estimation of the parton contributions to the nucleon spin. This may lead to a better 
understanding the spin components of the nucleon. 

In Refs. |16U17l[T^ll9j we investigated how those logarithmic corrections, in- 
fluenced the spin structure function, gi, at low values of Bjorken x. We formulated the 
evolution equations for the unintegrated parton distributions which included the complete 
resummation of the double logarithmic contributions, ln^(l/x). Afterwards, those equations 
were completed with the standard LO DGLAP evolution terms, in order to obtain the proper 
behaviour of the parton distributions at the moderate and large values of x. 

The equations obtained were applied to the following observables and processes: (i) to the 
nucleon structure function, gi, in the polarized deep inelastic scattering, (ii) to the structure 
function of the polarized photon, gj, in the scattering of a lepton on a polarized photon, and 
(iii) to the differential structure function, xj d'^gi/dxj dkj, in the polarized deep inelastic 
scattering accompanied by a forward jet. Case (iii) was proposed to be a test process for the 
presence and the magnitude of the ln^(l/x) contributions. For each process the consequences 
of including the logarithmic corrections were studied in a detail. 

Finally, some predictions for the observables, the asymmetry and the cross sections, in 
the processes (i)-(iii) were obtained. They were important to planned experiments with the 
polarized HERA and linear colliders, which would probe the region of low values of Bjorken 

X. 

2. Nucleon spin structure function, 

The total cross section for the polarized deep inelastic scattering is expressed by a contrac- 
tion of two tensors: the leptonic tensor, L^^, and the hadronic tensor, W^'^ , a ~ L^^W^^ PP. 
The hadronic tensor can be expanded in terms of the Lorentz invariants, and the expansion 
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coefficients define measurable structure functions ^Hj . The spin structure function, gi , is an 
expansion coefficient at the asymmetric spin term, e^^^" qpS^j^ where q denotes the momentum 
transfer, and 5 is a polarization vector of the nucleon. 

At low values of Bjorken x the asymptote of gi predicted by the Regge model shows a 
weak scaling with x |2Ull21j : 

5i(x,Q') = 7^(Q')x-"^^°\ (1) 

where is either a singlet {i = 5) or a non-singlet {i = NS) combination of the 

nucleon structure functions, g\ and g^. It is expected that as^Ns{^) ^ and that 05(0) 
"Afs(O) i-e. the singlet spin structure function with the Regge model is expected to show 
similar behaviour at low x as the non-singlet one [16j . 

The asymptotic form of gi obtained with the standard DGLAP evolution in QCD is more 
singular than this obtained with the Regge model jlUUllj . 

However, these predictions do not involve the contribution of the logarithms, ln(l/x), 
which is significant at low x. In the unpolarized DIS the resummation of the logarithms, 
ln(l/a;), was performed by the BFKL equation |22||2S|. The analysis of the polarized DIS 
showed that at low x the structure function of the polarized nucleon, 51 (x, Q^), was influenced 
by large double logarithmic corrections, Zn^(l/x), [TUllllj . Those corrections originated from 
the ladder diagrams and the nonladder (" bremsstrahlung" ) diagrams with the requirement 
of ordering the ratios, k'^/xn, for the exchanged partons |lUlllll[^l25j . The momentum, A;„, 
was the transverse momentum, and the fraction, x^, was the longitudinal momentum fraction 
of the nth parton exchanged. Ladder diagrams corresponded to the diagrams describing the 
forward scattering of a photon on a nucleon with a ladder of partons radiatively generated 
jlUlllllIT^] . Nonladder diagrams were the ladder diagrams of the forward scattering with 
some extra parton propagators attached. They were resummed with the infrared evolution 
equation pHI^ . 

Those logarithmic corrections were studied in a detail in Refs. jl6ll5j . Recursive equations 
for the resummation of the double logarithmic contributions were there formulated, both 
for the singlet and the nonsinglet component of gi. Those equations transformed to the 
moment space gave proper anomalous dimensions as derived from the infrared evolution 
equations mHillS]. 

The resummation of Zn^(l/x)was performed for the unintegrated parton distributions, 
fi {i = S,NS,g). Therefore the ordering requirement, k'^/xn < k'^^-^/xn+i, for the ex- 
changed partons was naturally included into the evolution equations. Moreover, the evolution 
/n^(l/x) was included into the equation kernels, and the nonperturbative input was separated 
as an inhomogeneous term. This allowed for including the kernels of the standard DGLAP 
evolution into the equations, which was necessary, in order to obtain a proper behaviour of 
the parton distributions at moderate and large values of x. 
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The full (unified) equation for the nonsinglet unintegrated parton distribution, /^rg was: 

( DGLAP ) 

, ~rr,2^^ f^dz fQy^dk\ fx . 



J X Z 



( Ladder 



/ qq 



1 dz fQy^ dk^ 
x z Jq2 k^ 

(Nonladder) 



(2) 



For a detailed form of the kernels see Appendix A. Matrices, Fg and Go, represent the octet 
partial waves and the colour factors respectively. They are described in a detail in Appendix 



A. Symbol 



(z) denotes the inverse Mellin transform of Fg/uj^ 



Fs/u;' 



f:Z-Fs{u)/u'. 

S—ioo ^TTl 



(3) 



with the integration contour located to the right of the singularities of the function Fg(Lj)/a;^. 

The singlet distribution, fs=\ ] , was resummed by a similar vector equation, formu- 

lated for the quark and the gluon components [TH] . 

The integrated parton distributions, Agj(x,Q^), were obtained with the unintegrated 
ones, using the relation: 



Ag,(x,Q2) = Agf (x) + J^^ ^ f,(^' = ^(1 + k'), 



(4) 



where the phase-space was extended from to W'^ = Q'^{l/x — 1) corresponding to the 
total energy squared measured in the center-of-mass frame ^[S]. The index, i = S, NS,g, and 
Agf (x) was the input parton distribution including the contributions from the nonpertur- 
bative region. 



Finally, the structure function, gi, at the LO accuracy was: 



9l{x, Q^) = ^ {AqNsix, Q^) + Aqsix, Q^)} 



(5) 



where (e^) = > ^^'^ denoted the number of active flavours {Nj = 3). 



The analysis of gi obtained with the unified evolution including the DGLAP terms at 
the LO accuracy was performed in jEl- The results obtained with the standard LO DGLAP 
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evolution and the unified evolution (DL+LO DGLAP) were similar at x > 10^'^. Signifi- 
cant discrepancies appeared at very low x, x < 10~^, when the corrections, l'n?{l/x) , were 
predominant (Fig. . 

For the gluon distribution, Ag, the discrepancies between the LO DGLAP and DL+LO 
DGLAP curves were large (Fig.l^J. This shows that the unified evolution of gluon distribution 
is driven by the ladder and the nonladder terms. 

Recently, the unified equation were completed by the NLO DGLAP terms |261l27j . The 
relations between the NLO DGLAP and the DL+NLO DGLAP curves of Refs. [2fill27| were 
similar as the relations between the LO DGLAP and the DL+NLO DGLAP curves obtained 



Unified evolution equations may be also applied to investigate the structure of the po- 
larized photon in a deep inelastic scattering of an electron on photon beams ^7] . The spin 
structure function of the polarized photon, gl, will be measured with future linear colliders 
e'^e~ and 67 [2312211301 • In particular, the scattering 67 will probe the gl at low values 
of Bjorken x. The partonic content of the polarized photon can be also measured with the 
process of the dijet photoproduction at the scattering of an electron on a proton j31j . 

Using the unified equations for the description of 67 scattering requires an extra inho- 
mogeneous term, A/cj(x,(5^), to be included into these equations. This term describes the 
pointlike coupling of the photon to quarks, antiquarks and gluons |28j . It joins the DGLAP 
evolution equations in the following way: 



where i = S, NS, g. There is no photon-gluon coupling at the LO of DGLAP, Akg{x, Q"^) = 

EHl. 



In Ref. |E1 we found an analytic solution of the simplified evolution equations. They 
described the DL evolution in an approximation, where the nonladder terms were neglected. 
We got the asymptotic form of gj at the limit, x ^ 0. As expected, gj scaled with x. The 
scaling coefficient was negative, and it was related in a simple way to the asymptotic scaling 
coefficients of the nucleon structure function, gi. 

Then we found a numerical solution for gl with the full evolution equations, taking into 
account the standard DGLAP evolution at the LO accuracy. The input parametrization 
fulfilled the sum rule for gj [H^IH^ : 



in dni. 



3. Spin structure of the polarized photon 



d 



Aqi{x,Q^) = Ak,{x,Q^) + as{Q^) {AP Aqi){x,Q^), 



(6) 



dlnQ^ 




(7) 
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We considered two limiting cases for the input parametrization of the quark and gluon dis- 
tributions: (i) the case, when Aq^^^ = and Ag^'^^ = 0, and the solutions of the evolution 
equations were generated radiatively from the inhomogeneous terms, Ak, (ii) the case, when 
Aq^^^ = but the gluon input was nonzero. In that case the input paraterization of gluon 
distribution was obtained with the vector- meson-dominance model (VMD), assuming the 
dominance of the mesons, p and w pill35j . 

Including the double logarithmic corrections significantly influcenced the behaviour of 
gj. The differences between the curves obtained with the LO DGLAP, NLO DGLAP and 
DL+LO DGLAP evolution strongly depended on the input parametrization of the parton 
distributions used (Figs. I3I4|) . It depended also on the infrared cut-off, (not shown). 
However, the estimated value of asymmetry, gJ/Fi was small, ~ 10"'^ at x = 10~^, i.e. it 
was difficult to measure. 

The gluon distribution, Ag, obtained strongly depended on the input parametrization 
used (Figs. I5I6|) . In case (i) Ag was negative at low values of x, in case (ii) it was positive. 
Future experiments will be then able to determine which parametrization is correct. 

4. Logarithmic corrections tested 

The presence and the magnitude of the logarithmic corrections, ln^{l/x) , can be tested 
in the polarized deep inelastic scattering accompanied by a forward jet. This idea was first 
applied for testing the BFKL resummation in the unpolarized DIS |361l37j . 

Assume, that a forward jet is produced in the DIS of an electron on a proton. The 
momentum coordinates of the jet are, {xj,kj), where xj is the longitudinal momentum 
fraction, and kj denotes the transverse momentum of the jet squared. The momentum 
transfer carried by the photon is q, and the Bjorken variable x is defined as: x = Q'^/{2pq), 
where = —q'^ and p is the momentum of the proton. 

If the jet coordinates fulfill the requirement: 

Xj » X, 

kj ~ (8) 

the jet is produced at low values of x/xj. A contribution to this process coming from 
the standard DGLAP evolution is then supressed {k'j ~ Q"^) [SHlEZlEHllMllinilll] , and the 
nonperturbative region is unpenetrated by the evolution, ln^(l/x) j42ll43j . 

In this process we measure the differential cross section. It is related to the momen- 
tum transfer, Q^, the electron energy fraction, y, and the differential structure function 
xjd'^gi/dxjdk'j. In analogy to the full structure function, gi, the differential structure func- 
tion may be expressed through the integrated parton distributions in a proton and the unin- 
tegrated distributions of quarks and antiquarks in a parton . 
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In Ref. we formulated the lv?{\/x) evolution equations for the unintegrated distribu- 
tions of (anti)quarks in a parton. In those equations we included the leading ladder terms, 
and neglected the subleading ones. We neglected also the DGLAP terms as the DGLAP 
evolution was supressed ^Sj- We found an analytic solution of those equations in the fol- 
lowing cases: (i) for the fixed coupling constant, as(/i^), where ^ = (fcj + Q'^)/2 |55lOU] . (ii) 
for the running coupling constant, 05 (/x^), where iJ? = fcj/C, and k'j^ denotes the transverse 
momentum of (anti)quark squared, and Q is its longitudinal momentum fraction. After nu- 
merical integrating the unintegrated distributions we got the differential structure function, 
xjd^gi/dxjdkj (Figs. [7|H|) . This function strongly changed with the value of the ratio, x/xj, 
and it depended also on the transverse momentum of the jet, kj. The strong dependence on 
x/xj was a direct effect of the logarithmic resummation. Comparing the results, we found 
that the case (ii) with the running coupling constant was more realistic than the case (i) 
with the fixed coupling constant. In both cases the Zn^(l/x) effects were much larger than 
the background described by the Born approximation (not shown). 

In Ref. jl9j we estimated the cross section and the asymmetry for the DIS accompanied 
by a forward jet. We used the kinematical cuts as planned for the future collider, the polar- 
ized HERA j44ll45lll^ . The effects of /n^(l/a;) resummation contributed significantly to the 
asymmetric cross section da/dx, where x was the Bjorken variable (Fig. |^. The asymmetry 
obtained was, however, small, and changed between —0.01 and —0.04 at low x. 

5. Summary 

We performed the analysis of the polarized deep inelastic scattering in the region of 
low values of Bjorken x. We used the evolution equations formulated for the unintegrated 
parton distributions. Those equations included the complete resummation of the double 
logarithmic contributions, ln^(l/x), and the standard LO DGLAP evolution terms. The 
DGLAP evolution was included into those equations, in order to obtain the proper behaviour 
of the parton distributions at the moderate and large values of x. 

Using that formalism, we obtained predictions for the following observables: (i) the nu- 
cleon structure function, gi, in the polarized deep inelastic scattering, (ii) the structure 
function of the polarized photon, gj, in the scattering of a lepton on a polarized photon, 
and (iii) the differential structure function, xj d'^gi/dxj dkj, in the polarized deep inelastic 
scattering accompanied by a forward jet. 

For each process the consequences of including the logarithmic corrections were studied in 
a detail. Generally, we observed an enhancement of the magnitude of g'l, gj and Ag obtained 
with the unified evolution as compared to the pure LO DGLAP results. Especially, the be- 
haviour of gluon distributions with the unified evolution at low x was clearly dominated by 
the large contribution of the ladder and nonladder terms. Some predictions for the observ- 
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ables, the asymmetry and the cross sections, in the processes (i)-(iii) were also obtained. The 
case (iii) tested the presence and the magnitude of the ln^(l/a;) contributions in DIS. In that 
case the asymmetry parameter estimated with the kinematical cuts of the future cohider, 
polarized HERA, was between —0.04 and —0.01 at low x. 

To sum up, our observations suggest that the standard DGLAP evolution is not complete 
at low X, where the effects of ln'^{x) resummation are large and therefore cannot be neglected. 
This result is important to planned experiments with the polarized HERA and linear colliders, 
which would probe the region of low values of Bjorken x. 
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Appendix A 

Here a brief description of the evolution kernels of Eq. @ is given. DGLAP kernels were 
taken from Ref. pp. Here the DGLAP kernel AP includes only the LO terms: 

AP = AP(°). (9) 

The ladder kernels corresponding to the LO DGLAP kernels at the longitudinal mo- 
mentum transfer, z = 0, ^H] generate the double logarithmic corrections in the region of 
<k^< Q^/z. 

The nonladder kernels were obtained in Ref. ^5] from the infrared evolution equations 
written for the singlet partial waves Fq, Fg |10lllllOHl25j . In we noticed that extending 
the kernel of the double logarithmic evolution equations from the ladder one, 

as{Q^)APgg/Lo, (10) 

to the modified one, 

STsiQ^) [aPJu; - (F8(c^) Go^/ {27:^oj^)) , (11) 
gave a proper anomalous dimension as derived from the infrared evolution equations. 
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Matrix Gq contained colour factors resulting from attaching the soft gluon to external 
legs of the scattering amplitude : 



Ar2-1 



Gn 







2N 

N 



(12) 



where was the number of colours. 

Further, it was checked that the Born approximation of Fg, 

F3^-"M«8vr25;(Q2)^. 



(13) 



gave accurate results for the DL evolution. Martix Mg was a splitting function matrix in the 
colour octet channel, 



/ L -EiL 

( 2N 2 

\ N 2N 
The inverse Melin transform of Fq"^"'{u!) then read : 

-pBorn 



(14) 



{z) =4TT'asiQ')M8ln'{z). 



(15) 



The evolution equation Q includes the nonladder corrections in the Born approximation 
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Figure 1. Structure function, gi{x,Q^), at = 10 GeV^ plotted as a function of x. Solid line corresponds 
to the results with the full ln'^(l/x) resummation, where the nonladder terms, ladder terms and the DGLAP 
kernels were included. Dotted line shows the pure DGLAP evolution. Thin solid line shows the nonperturbative 
input, pS^*"'"' , and dashed line shows the incomplete DL resummation, where the nonladder terms were neglected 
and the DGLAP terms were included. 
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Figure 2. The spin dependent gluon distribution, Ag{x,Q^), at = 10 GeV^ plotted as a function of x. 
Solid line corresponds to the results with the full ln'^{l/x) resummation, where the nonladder terms, ladder 
terms and the DGLAP kernels were included. Dotted line shows the pure DGLAP evolution. Thin solid line 
shows the nonperturbative input, Ag^°\ and dashed line shows the incomplete DL resummation, where the 
nonladder terms were neglected and the DGLAP terms were included. 
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Figures. Structure function, gJ{x,Q'^)/aem, at = 10 GcV'^ obtained after solving the unified evolution 
equations with the input pararnctrization (i). Solid line corresponds to the results with the full ln^(l/x) 
resurnrnation, where the nonladder, ladder corrections and the LO DGLAP kernels were included. Dashed 
line shows the LO DGLAP evolution, and dotted line shows the NLO DGLAP evolution. 
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Figure4. Structure function, gJ{x,Q'^)/aem, at = 10 GeV^ obtained after solving the unified evolution 
equations with the input parametrization (ii). Solid line corresponds to the results with the full ln'^{l/x) 
resummation, where the nonladder, ladder corrections and the LO DGLAP kernels were included. Dashed 
line shows the LO DGLAP evolution, and dotted line shows the NLO DGLAP evolution. 
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Figures. Spin dependent gluon distribution in photon, IS.g'^ {x,Q^) / cteri-i . at = 10 GeV'^ obtained after 
solving the unified evolution equations with the input paranietrization (i). Sohd hue corresponds to the results 
obtained with the full ln^{l/x) rcsuniniation, where the nonladder, ladder corrections and the LO DGLAP 
kernels were included, dashed line shows the LO DGLAP evolution. 
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Figure 6. Spin dependent gluon distribution in photon, Ag'^{x,Q'^)/aem, at = 10 GeV^ obtained after 
solving the unified evolution equations with the input parametrization (ii) . Solid line corresponds to the results 
obtained with the full ?n^(l/a;) resummation, where the nonladder, ladder corrections and the LO DGLAP 
kernels were included, dashed line shows the LO DGLAP evolution. 
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Figure?. The differential spin structure function, xj q^^qj^^ , for the fixed Qs (case (i)) plotted as a function 
of the longitudinal momentum fraction carried by a jet, xj. We show predictions for the three different values 
of the transverse momentum of the jet squared, kj = 5 GeV'^, 10 GeV"^ and 20 GeV^. Those calculations were 
performed at = 10 GeV^ and x = 10"*. 
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Figures. The differential spin structure function, xj g^^gj^i , for the running cis (case (ii)) plotted as the 
function of the longitudinal momentum fraction carried by a jet, xj. We show predictions for the three 
different values of the transverse momentum of the jet squared, kj =5 GeV^, 10 GeV^ and 20 GeV^. Those 
calculations were performed at = 10 GeV^ and x = 10~*. 
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Figure 9. The cross-section, for the forward jet production in the polarized deep inelastic scattering. 
Figures 2a and 2b show the results with the double logarithmic in'^(l/^) effects included. They correspond 
to two choices of the scale i?: (i) = {kj + Q'^)/2 (Fig. 2a), (ii) /i^ = k}/^ (Fig. 2b). Fig. 2c shows the 
cross-section ^ in the Born approximation where the double logarithmic resummation was neglected. 



